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ABSTRACT 

We study the quasi-simultaneous near-IR, optical, UV, and X-ray photometry 
of eleven 7-ray selected blazars for which redshift estimates larger than 1.2 have 
been recently provided. Four of these objects turn out to be high-power blazars 
with the peak of their synchrotron emission between ~ 3 x 10"'^^ and ^ 10^^ 
Hz, and therefore of a kind predicted to exist but never seen before. This dis- 
covery has important implications for our understanding of physical processes 
in blazars, including the so-called "blazar sequence", and might also help con- 
straining the extragalactic background light through 7-ray absorption since two 
sources are strongly detected even in the 10—100 GeV Fermi -LAT band. Based 
on our previous work and their high powers, these sources are very likely high- 
redshift flat-spectrum radio quasars with their emission lines swamped by the 
non-thermal continuum. 

Key words: BL Lacertac objects: general — radiation mechanisms: non- 
thermal — gamma-rays: observations — radio continuum: galaxies — X-rays: 
galaxies 
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1 INTRODUCTION 

Blazars are radio loud active galactic nuclei (AGN) with 
their jets pointing towards the observer (see e.g. Bland- 
ford & Rees 1978; Urry & Padovani 1995). Observation- 
ally, they come in two main flavours with very different 
optical spectra. Namely, Flat Spectrum Radio Quasars 
(FSRQs) have strong, broad emission lines, just like ra- 
dio quiet quasars, while BL Lacertae objects (BL Lacs) 
show at most weak emission lines, sometimes display ab- 
sorption features, and can also be completely featureless. 

Blazars emit variable, non-thermal radiation across 
the whole electromagnetic spectrum, which includes two 
components forming two broad humps in a. v — vfv rep- 
resentation. The low-energy one is attributed to syn- 
chrotron radiation and the high-energy one is usually 
thought to be due to inverse Compton radiation (e.g. 
Abdo et al. 2010b). The peak of the synchrotron emis- 
sion (I'poak) can occur over a wide range of frequencies, 
from about ~ 3 x 10^^ Hz to over 10^* Hz, reflecting 
the maximum energy at which particles can be acceler- 
ated (e.g. Giommi et al. 2012b). Blazars with rest-frame 
^'pcak < 10^^ Hz are called Low Synchrotron Peaked 
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(LSP) sources, while those with 10^* Hz < 



peak 



lO'"" Hz and v, 



peak 



10^^ Hz are called Intermediate 
and High Synchrotron Peaked (ISP and HSP) sources, re- 
spectively (Abdo et al. 2010b). These deflnitions expand 
on the original division of BL Lacs into LBL and HBL 
sources first introduced by Padovani & Giommi (1995). 

The location of the synchrotron peak implies dif- 
ferent mechanisms for the X-ray emission, namely in- 
verse Compton in LSPs, which have hard energy index 
(ctx ^ 1) spectra, and an extension of the synchrotron 
emission responsible for the lower energy continuum in 
HSPs, which typically display steeper (qx ~ 1.5) X-ray 
spectra^ (e.g. Padovani & Giommi 1996; Wolter et al. 
1998; Padovani et al. 2004). 

The distribution of the synchrotron peak frequency 
is extremely different for the two blazar classes. While 
the rest-frame i-'p^ak distribution of FSRQs is strongly 
peaked at low frequencies {(Vpcak) ~ lO^'^'^*"'^ Hz) and 
never reaches high values (Vpazk ~ 10^*'^ Hz) (Giommi et 
al. 2012b), the 

^peak distribution of BL Lacs is shifted to 



^ Unless I'pjjg^jj is in the X-ray band, where one samples the 
top of the synchrotron hump and the X-ray spectrum is also 
relatively flat (Padovani & Giommi 1996). 
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higher values by at least one order of magnitude and can 
reach i^p^ak ^ high as ;i 10^* Hz. In other words, while 
BL Lacs can be of the LSP and HSP type, FSRQs appear 
to be solely LSPs. 

The origin of this difference has been discussed at 
length in the literature. In a series of papers, Padovani, 
Giommi, and collaborators, looked for HSP FSRQs using 
multi-wavelength catalogues and the deep X-ray radio 
blazar survey (DXRBS; e.g. Padovani et al. 2003) with 
only limited success. In fact, while FSRQs with Vp^^y^ pos- 
sibly larger than previously known were discovered, 
follow-up X-ray observations showed in most cases a rela- 
tively flat X-ray spectrum. This, together with their over- 
all spectral energy distributions (SEDs), was suggestive 
of inverse Compton emission (e.g. Padovani et al. 2002; 
Landt et al. 2008) and did not support a possible HSP 
classification. The only exception is RGB J1629-I-4008, 
with t'peak ~ 2 X 10^'' Hz and steep X-ray spectrum 
(Padovani et al. 2002). ROXA J081009.9-H384757.047 
was also suggested to be an HSP (Giommi et al. 2007) 
but new WISE data (Wright et al. 2010) indicate that 
the near-IR and optical flux can be attributed to accre- 
tion (the big "blue bump") and that the peak of the 
non-thermal synchrotron radiation is located in the far 
infrared, making this object a typical LSP. 

The purported non-existence of HSP FSRQs implied 
also the lack of high-power HSP blazars (since FSRQs are 
more powerful than BL Lacs), which was also hinted at 
by a suggested anti-correlation between bolometric lumi- 
nosity and !-'pcak (the so-called "blazar sequence"; Fos- 
sati et al. 1998). This has been interpreted by some re- 
searchers as an intrinsic physical difference between LSPs 
and HSPs. Since HSPs are typically characterized by a 
low intrinsic power and external radiation fleld, given 
their very weak or absent emission lines, cooling was 
thought to be less dramatic than in LSPs allowing parti- 
cles to reach energies high enough to produce synchrotron 
emission well into the X-ray band (Ghisellini et al. 1998). 
In this respect, RGB 1629-1-4008, with its relatively low 
radio power, Psghz ~ 4 x 10^* W Hz~^, was regarded as 
fitting on the sequence as well (Padovani et al. 2002). 

Giommi et al. (2012a) have recently presented a new 
scenario, where blazars are classifled as FSRQs, BL Lacs, 
LSPs, or HSPs according to a varying mix of the Doppler 
boosted radiation from the jet, emission from the accre- 
tion disk, the broad line region, and the host galaxy. This 
hypothesis predicts the existence of high power - high 
^^poak blazars, which however would be very hard to iden- 
tify because of their featureless optical spectra and, there- 
fore, lack of redshift. This is because when both fp^ak a-nd 
radio power are large the dilution by the non-thermal 
continuum becomes extreme and all optical features are 
washed away. 

Very recent results provide a simple way to test 
this prediction and look for the existence of previously 
unidentified high power HSPs. Rau et al. (2012) have pre- 
sented redshift constraints for 103 blazars from the Fermi 
2LAC catalogue (Ackermann et al. 2011) by fitting SED 
templates to their UV-to-near-IR multi-band photome- 
try obtained quasi-simultaneously with Swift/WOT and 
GROND (Greiner et al. 2008). In fact, attenuation due 
to neutral hydrogen along the line of sight at the Lyman 
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Figure 1. The near-IR through UV SEDs for the eleven 
blazars with Zp^ot > 1-2 from Rau et al. (2012) in a u — ufu for- 
mat. The lines are the best fit SEDs with red, thick lines indi- 
cating our HSPs. From bottom to top: CRATES J0402-2615, 
CRATES J1312-2156, OM 235, CRATES J0630-2406, 
SUMSS J053748-571828, PKS 0600-749, PKS 0332-403, 
CLASS J2352-I-1749, PKS0047+023, PKS 0055-328, RX 
J0035. 2-1-1515. 

limit allows an accurate estimate of the redshift of the ab- 
sorber and thereby provides a reliable lower limit to the 
blazar redshift. Eleven of their objects have Zphot > 1-2 
and for eight such sources these are the first reliable red- 
shift measurements. Since the spectrum of most of these 
sources is featureless and high redshift typically means 
high power, this is an ideal sample to look for high power 
HSPs with their lines completely swamped by the non- 
thermal continuum, which is the purpose of this work. 

Throughout this paper we use a ACDM cosmology 
with Ho = 70 km s"^ Mpc"\ Qm = 0.27 and Ha = 0.73 
(Komatsu et al. 2011). 



2 SEDS OF THE NEWLY IDENTIFIED 
HIGH REDSHIFT BLAZARS 

Fig. 1 plots the near-IR through UV SEDs of all eleven 
blazars with Zphot > 1.2 from Rau et al. (2012) in a 
u — !^/^ representation. Taking into account the obvi- 
ous decrease in fiux at high frequencies due to neu- 
tral hydrogen absorption, it is already clear from these 
SEDs covering only about one decade in frequency that 
a few of the sources do appear to have i^p^ak ~ 10^^ Hz, 
which qualifies them as HSPs. We have then assembled 
the SEDs of all eleven sources over the whole electro- 
magnetic spectrum by using the SED builder tool avail- 
able at the ASDC (tools.asdc.asi.it/SED, Stratta et al. 
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Figure 2. The rest-frame SED of RX J0035.2+1515. Red 
points represent the quasi-simultaneous GROND and Swift 
data. The hght red 7-ray points are from the second Fermi 
-LAT catalog (Abdo et al. 2012). Gray points are non- 
simultaneous archival data. The solid line is a third degree 
polynomial fit to the radio points, simultaneous GROND and 
Swift UVOT data not afl'ected by UV absorption, and XRT 
data. 



2011). As mentioned in the Introduction, a key require- 
ment for an HSP classification is a steep X-ray spectrum 
(for I'peak ~ 10^ Hz). Since the Swift/X-r&y Telescope 
(XRT) observes contemporaneously with Swift/WOT , 
our sources have quasi-simultaneous coverage from the 
GROND bands into the X-ray regime, i.e. between 10^* 
and 10^® Hz. This is an ideal situation to determine 
for HSPs. 

It turns out that four sources are HSPs, whose main 
observational parameters are given in Table 1. We deter- 
mined rest-frame i^p^ak values by fitting a third degree 
polynomial^ to the radio data, the GROND and UVOT 
data below the Lyman limit, and the quasi-simultaneous 
Swift/XRT data. We also fitted the X-ray data (using the 
XSPEC12.0 package) adopting a simple power law spec- 
tral model with the low energy absorption fixed to the 
Galactic value in the direction of the source (Dickey & 
Lockman 1990). In two cases, given the relatively short 
exposure times, we used the XIMAGE package (Giommi 
et al. 1992) V4.5^ to estimate the source count rates in 
the soft (0.5 - 2.0 keV) and hard (2.5-10 keV) X-ray en- 
ergy bands and derive an estimate of the spectral slope. 
In the case of CRATES J1312-2156, however, we could 
only get a 3a upper limit at higher energies, which trans- 
lates into a lower limit to the slope. 

The SEDs for the HSPs are shown in Figs. 2-5 
while Table 2 gives z^pcak, -f-pcak, A.4GHz, ctx, and a^. All 
four sources are characterized by t'pcak > 3 x 10^^ Hz and 
a rising 7-ray spectrum {a-, ^ 1). Three of the sources 
display also a very steep X-ray spectrum (qx ^ 1-6), 
while for the fourth one we can only say that Ox > 0.5. 
ipcak and Pi.4GHz are quite high, being > 4 x lO**® erg 
s-i and > 2 X lO^"^ W Hz'S i.e. typical of FSRQs. We 



^ This represents the simplest way to estimate i^^^^y^ in a 
straightforward and model-independent way. 
3 XIMAGE is part of the HEASARC Xanadu standard soft- 
ware package for multi-mission X-ray astronomy. 
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Figure 3. The rest-frame SED of SUMSS J053748- 

Data points and solid line are as in Fig. 2. 



571828. 
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Figure 4. The rest-frame SED of CRATES J0630-2406. Data 

points and solid line arc as in Fig. 2. 



also note that SUMSS J053748-571828 shows clear evi- 
dence of X-ray variability of the type displayed by HSPs 
in the X-ray band when this is beyond t'peak- Similarly, 
CRATES J0630-2406 also shows evidence of X-ray vari- 
ability, with Swift/X.KT data taken at different times still 
displaying a steep X-ray spectrum. 
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Figure 5. The rest-frame SED of CRATES J1312-2156. Data 
points and solid line are as in Fig. 2. 
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Table 1. Observational parameters 



Name Fermi Name Optical Coordinates'^ ^phot 

/l.4GHz f"^ /2-lOfceV /l-lOOGeV 

mjy AB erg/cm-^ /s 10"-'^'^ ph/cm^/s 



RX J0035. 2+1515 


2FGLJ0035. 2+1515 


00 35 14.7 


+ 15 15 04.2 


1 90+0.14 


18.8 


17.10 


3.3 10- 


ts 


10.8 


SUMSS J053748-571828 


2FGLJ0537.7-5716 


05 37 49.0 


-57 18 30.4 


1 cc;+0.09 

^■'^^-o.is 


99.8^ 


17.48 


5.5 10- 


13 


3.9 


CRATES J0630-2406 


2FGLJ0630.9-2406 


06 30 59.5 


-24 06 46.2 


1 60+°" 
^■°^-0.05 


105.8 


15.90 


10.7 10 


-13 


33.8 


CRATES J1312-2156 


2FGLJ1312.4-2157 


13 12 31.6 


-21 56 23.5 


1 7Y+0.09 
-^■"-0.11 


349.8 


16.56 


5.1 lo- 


13 


21.7 



''From Rau et al. (2012) 
'Radio flux density at 843 MHz 



Table 2. Inferred parameters 



Name 


log , 

^ peak 

Hz (rest-frame) 


log .J^pcak 

erg s--"- 


log Pl.4GHz 
W Hz^i 


fix 

(0.5 - 10 kcV) 


(100 MeV - 100 GcV) 


RX J0035. 2+1515 
SUMSS J053748-571828 
CRATES J0630-2406 
CRATES J1312-2156 


15.8 
15.5 
16.0 
15.6 


46.6 
46.6 
47.1 
46.9 


26.2 

27.1'' 

27.1 

27.7 


1 9+"-i 
1 &+°-^'' 

^^■"-0.3 
1 9+0.2 

> 0.5^= 


0.61 + 0.11 
0.73 + 0.21 
0.79 + 0.06 
1.02 + 0.07 



"From Abdo et al. (2012) 

''Extrapolated from 843 MHz assuming Or = 

•^Derived from the count rates in the 0.5 — 2.0 and 2.5 — 10 keV bands 
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Figure 6. The bolometric power (represented by the sum of 
the synchrotron and inverse Compton peak powers) of the flux- 
limited samples of blazars of Giommi et al. (2012b) is plotted 
against i^p^j^j. (see Giommi et al. 2012b, for details). The four 
objects discovered in this paper, represented by the red filled 
points, are clearly located in an area of the diagram that was 
not populated in the past. 



3 DISCUSSION AND IMPLICATIONS 

All of the observational evidence suggests that we have 
discovered four blazars of a previously unknown type, 
that is sources having both high power (Lpcak ~ 6 x 
lO''^ erg s"^) and high i/^^^^ (« 5 x 10^^ Hz). These 
objects have also radio powers 2 — 3 orders of magnitude 
larger than the only previously known HSP FSRQ and 
therefore are firmly into the FSRQ regime (unlike RGB 
1629+4008, which was near the low-luminosity end of the 
FSRQ luminosity function: Padovani et al. 2002). 

To put things into perspective. Fig. 6, adapted from 
Giommi et al. (2012b), plots the bolometric power (i.e.. 



the sum of the synchrotron and inverse Compton peak 
powers) for blazars selected in the radio. X-ray and 7- 
ray-bands, against i^peak- Our HSPs end up in an empty 
region of the plot, where no objects were found before. 
Similar results are obtained if all 2LAC BL Lacs with 
redshifts are also included. 

Our i^peak estimates are robust, as they are based 
on quasi-simultaneous coverage over four decades in fre- 
quency and are corroborated in at least three cases by a 
steep X-ray spectrum. Contamination by a thermal (ac- 
cretion disk) component is ruled out by the featureless 
spectra of these sources and also by the fact that the flux 
of an accretion disk fitted using the scaling of Giommi et 
al. (2012a) would be well below the observational data. 
Our sources also end up in the HBL region of the - 
!^pcak plot of Padovani et al. (2004) (their Fig. 5) . 

These four sources are also the only ones with a~, 5C 
1, as the other seven in the Rau et al. (2012) sample have 
> 1 (Abdo et al. 2012). While one is borderline with 
~ 1, three have a~i 0.79. We note that no FSRQ, 
LSP, or ISP reach these low values in the Clean Sample 
of Ackermann et al. (2011). It is also interesting to note 
that 4/6 of the remaining objects in the Rau et al. (2012) 
sample (excluding the FSRQ OM 235) have v^^^^ > lO" 
Hz and therefore qualify as ISPs. 

In the radio power - i^p^ak plot of Giommi et al. 
(2012a) (their Fig. 11) the newly discovered HSPs end 
up in the region of parameter space mostly populated 
by sources without redshift, which agrees with the fact 
that, indeed, these sources have totally featureless optical 
spectra. Based on that paper, our sources are very likely 
FSRQs with their emission lines swamped by the non- 
thermal continuum. Given their large powers, this makes 
sense also on statistical grounds. Indeed, for Lboi > 5 x 
10*'' erg s-^, which is the smallest value for our sources, 
previously known FSRQs outnumber BL Lacs by a factor 
> 10 in Fig. 6. 



© 2012 RAS, MNRAS 000, 1-5 



High power - high synchrotron peak blazars 5 



The discovery of high power HSPs has important 
imphcations and consequences for a variety of issues. We 
mention here some of the most prominent ones: 

(i) our HSPs end up in an empty region of the Lpcak - 
''peak plot of Meyer et al. (2011) (see their Fig. 4). In fact, 
the closest sources with t'peak similar to those of our ob- 
jects have Lpcak more than an order of magnitude smaller. 
This implies that the so-called "blazar sequence" (Fossati 
et al. 1998; Ghisellini et al. 1998) is indeed heavily influ- 
enced by the selection effects described by Giommi et al. 
(2012a); 

(ii) it is well established that a large fraction of BL 
Lacs (- 43% in BZCAT Massaro et al. (2009) and 
> 50-60% of the BL Lacs in the Fermi 1 and 2 year AGN 
catalogs: Abdo et al. 2010a; Ackermann et al. 2011) have 
no measured redshift due to the lack of any detectable 
feature in their optical spectrum, despite the use of 8/10- 
m class optical telescopes for the spectroscopy identifica- 
tion campaign. This paper, together with the results of 
Rau et al. (2012), confirms the suggestion by Giommi et 
al. (2012a) that these sources are high-redshift blazars, 
most likely FSRQs with their emission lines swamped by 
the jet. Therefore: a) one should not assume for these ob- 
jects a redshift typical of other BL Lacs but a much larger 
one (typically 2; > 1); b) these sources should be included 
with the FSRQs when studying cosmological evolution 
and luminosity functions, since their exclusion is bound 
to bias the results; 

(iii) the study of distant very high energy (VHE) emit- 
ting AGN is of fundamental importance to constrain the 
extragalactic background light (EBL). The 7-ray flux 
from distant blazars is in fact absorbed on its way from 
the source through its interaction with EBL photons. One 
can then use the measurement of the induced distortions 
to derive constraints on the EBL density (e.g. Mazin & 
Raue 2007). The discovery of high redshift, 7-ray emit- 
ting blazars, strongly detected in the Fermi highest en- 
ergy band (10- 100 GeV) in the case of RX J0035.2-M515 
and CRATES J0630-2406 (see Figs. 2 and 4), is there- 
fore very important also in this context. Some of the 
high-redshift sources of the type described in this pa- 
per may be detectable by existing Cherenkov telescopes, 
as the 10 — 100 GeV Fermi fluxes of both sources dis- 
cussed above are within the range of observed values of 
blazars that have already been detected at TeV energies. 
The next generation of VHE telescopes, particularly the 
Cherenkov Telescope Array, will be able to detect many 
such sources, especially if the low energy threshold will 
reach ~ 10 GeV, as planned. 
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